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The full gluon propagator relevant for the description of the truly non-perturbative QCD dynam- 
ics, the so-called intrinsically non-perturbative gluon propagator has been derived in our previous 
work. It explicitly depends on the regularized mass gap, which dominates its structure at small 
gluon momentum. It is automatically transversal in a gauge invariant way. It is characterized by 
the presence of severe infrared singularities at small gluon momentum, so the gluons remain mass- 
less, and this does no depend on the gauge choice. In this paper we have shown how precisely the 
renormalization program for the regularized mass gap should be performed. We have also shown 
how precisely severe infrared singularities should be correctly treated. This allowed to analytically 
I formulate the exact and gauge- invariant criteria of gluon and quark confinement. After the renor- 

f~* > ) , malization program is completed, one can derive the gluon propagator applicable for the calculation 

of physical observables, processes, etc., in low-energy QCD from first principles. 
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I. INTRODUCTION 



O i 1 In our previous work Jj the mass gap responsible for the non-perturbative dynamics in QCD has been introduced 
JL \ through the definition of the subtracted full gluon self-energy. So it is defined as the value of the regularized full 
gluon self-energy at some finite point. The mass gap is mainly generated by the nonlinear interaction of massless 
gluon modes. We have explicitly shown that QCD is the color gauge invariant theory at non-zero mass gap as well. 
All this allows one to establish the structure of the full gluon propagator in the explicit presence of the mass gap 
In this case, the two independent general types of the formal solutions for the regularized full gluon propagator 
have been found (2j. No truncations/ approximations/assumptions, as well as no special gauge choice are made for the 
^vq ' regularized skeleton loop integrals, contributing to the full gluon self-energy. The so-called massive solution, which 
leads to an effective gluon mass has been established. The general nonlinear iteration solution for the full gluon 
propagator depends explicitly on the mass gap. It is always severely singular in the q 2 — ► limit, so the gluons 
remain massless, and this does not depend on the gauge choice. However, we have argued that only its intrinsically 
[~^. ■ non-perturbative (INP) part, defined by the subtraction of all the types of the perturbative theory (PT) contributions 
("contaminations") from the general nonlinear iteration solution, is interesting for confinement. In this way it becomes 
0^ . automatically transversal in a gauge invariant way and excludes the free gluons from the theory. This is important 
for correct understanding of the color confinement mechanism, as underlined in our previous works [l], 0] • 

The INP gluon propagator is to be used for the numerical calculations of physical observables, processes, etc. in 
low-energy QCD from first principles. But before the two important problems remain to solve. The first problem is 
how to perform the renormalization program for the regularized mass gap A 2 = A 2 (A,ct,£, g 2 ), and to see whether 
the mass gap survives it or not. The second problem is how to treat correctly severe infrared (IR) singularities 
(q 2 )~ 2 ~ k , k = 0,1,2,3,... inevitably present in this solution (let us remind that the signature is Euclidean). Just 
these problems are to be addressed and solved in this paper. 

II. THE INP GLUON PROPAGATOR 

The full gluon propagator, which is relevant for the description of the truly NP QCD dynamics, the so-called INP 
gluon propagator derived in Ref. [1] is as follows: 

D™ P (q) = iT, v {q)d INP {q 2 )^ = i T, v (q)^L(q 2 ), (2.1) 



* gogohia@rmki.kfki.hu 



2 

and 

oo -2 00 » 2 00 

L( g ») S X( fl »;A») = X;(|-) ** = E(|") E«W ( 2 - 2 ) 

fc=0 /c=0 m=0 

where A 2 = A 2 (A, a,£,g 2 ) and the residues $fc = $fc(A, a,£,g 2 ) and thus $fcm = <£fc m (A, a, £, g 2 ) as well. 

Let us recall some interesting features of the INP gluon propagator (2.1) 2]. First of all, it depends only on the 
transversal degrees of freedom of gauge bosons. Also, its functional dependence is uniquely fixed up to the expressions 
for the residues and it is valid in the whole energy/momentum range. It explicitly depends on the mass gap, so 
that when it formally goes to zero (the formal PT A 2 — limit) this solution vanishes. Also, it is characterized by 
the explicit presence of severe (i.e., NP) IR singularities (q 2 )~ 2 ~ k , k = 0,1,2,3,... only. The sum over m indicates 
that an infinite number of iterations (all iterations) invokes each severe IR singularity labeled by k. Apart from the 
structure (A 2 /{q 2 ) 2 it is nothing but the corresponding Laurent expansion. This solution is completely free of all 
the types of the PT contributions ("contaminations"), and thus is exactly and uniquely separated from the PT gluon 
propagator, indeed. Due to the regular dependence on the mass gap, it dominates over the free gluon propagator in 
the deep IR (q 2 — » 0) limit and it is suppressed in the deep ultraviolet (UV) (q 2 — > 00) limit (let us remind that the 
free gluon propagator behaves like 1/q 2 in the whole momentum range). In the rest of this paper we will omit the 
superscript "INP", for simplicity. We will restore it if it will be necessary for clarity of our discussion. 



III. THE GENERAL MULTIPLICATIVE RENORMALIZATION PROGRAM 

One of the remarkable features of the solution (2.1)-(2.2) is that its asymptotic at infinity (q 2 — > 00) is to be 
determined by its A 2 /(q 2 ) 2 structure only, since all other terms in its expansion are suppressed in this limit. It is well 
known that such a behavior at infinity is not dangerous for the renormalizability of QCD [H, [H, 0, [(| . However, the 
regularized mass gap itself is of the NP origin. So its renormalization in order to get finite result in the A — > 00 limit 
cannot be done in standard way. Another main problem, closely connected to the mass gap, is its violent structure in 
the deep IR region (q 2 — ► 0). Fortunately, there exist two mathematical theories which are of a great help in this case. 
The distribution theory (DT) 0], which should be complemented by the dimensional regularization method (DRM) 
Q (some important issues of these theories are present and discussed in appendix A in more detail). 

For further purpose, let us present our solution (2.1) as follows: 



D^(g)=iT MV (g)V(AT +i (<n 



. 2\k+l, 
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(3.1) 



As it has been emphasized in appendix A, a special regularization expansion is to be used in order to deal with the 
severe IR singularities (q 2 )~ 2 ~ k , since standard methods fail to control them. If q 2 is an independent loop variable, 
then the corresponding dimensionally regularized expansion is given in Eq. (A10), namely 
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(3.2) 



Substituting it into the previous expansion, one obtains 



D^{q) = iTM- E( A2 ) fc+1 *fc \a(k)[S 4 (q)] {k) + O k (e) 



k=0 



(3.3) 



so instead of the Laurent expansion (2.1) we obtain dimensionally regularized (in powers of e) expansion for the INP 
gluon propagator. In the presence of the IR regularization parameter e all the quantities and parameters depend, 
in general, on it. For further purpose it is convenient to introduce the following short-hand notations, namely 
A 2 = A 2 (A, e) and = 3>fc(A, e) and omit the dependence on the other parameters, for simplicity. Evidently, the 
renormalization of the mass gap only is what that matters, since $fe(A, e) is the dimensionless quantity (see below). 
So let us define the UV and IR renormalized mass gap as follows: 



A 2 (X;e) = Z(X;e)A 2 R , 



00, 



(3.4) 
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where Z = Z(\; e) is the multiplicative renormalization (MR) constant of the mass gap. Here and everywhere below 
all the quantities with sub- or superscript "R" denote the quantities which exist in the A — > oo limit, i.e., by definition, 
they are UV renormalized (and their dependence on other parameters is to be neglected, apart from the coupling 
constant, which can be "running" itself, see section V below). At the same time, all the quantities with bar are IR 
renormalized, i.e., by definition, they exist as e — > + . That is why here and everywhere below in all the renormalized 
quantities the dependence on A and e is not shown explicitly, as well as on other possible parameters. The factorization 
of the corresponding MR constants in order to get the finite results separately in the A — ► oo and e — > 0+ limits is 
a particular case of the general MR program in order to remove the dependence on A and e in their corresponding 
limits. The NP (quadratic) UV divergences present in the regularized mass gap are in close intrinsic link with the 
NP IR singularities, which both are present in the INP QCD gluon propagator (2.1). 
Substituting the relation (3.4) into the previous expansion (3.3), one arrives at 



D^(q) = iT^(q)- Y^(A 2 R ) k + 1 Z k + 1 (\; e)$ fc (A; e) \a(k)[6 4 (q)]^ + O k (e) 



0+. 



(3.5) 



k=0 



The only way to remove the pole 1/e from the previous expansion and get the final results after the removal the 
dependence on A and e is to put 



Z fe+1 (A;e)$ fc (A;e) = efi£, e^0+, A ^ oo, fc = 0,l,2,3, 



(3.6) 



since all the NP IR singularities labeled by k are independent from each other. This is the general convergence 
condition which makes the gluon propagator free of A and e in their corresponding limits. 
So the previous expansion becomes 



D«M =iT^q)^(Al) k+1 nx\a(k)[5\q)]W + Q k (e) 



k=0 



0+ 



(3.7) 



Through the chain of the relations fig = ] fe+1 , Af^f = (A^) fe , one obtains 



0+ 



(3.8) 



which shows that one can start from the fc-dependent regularized mass gap in Eq. (3.3), but, nevertheless, coming to 
the same final Eq. (3.8), containing the fc-dependent renormalized mass gap. 

The renormalization of the mass gap itself has been defined by the relation (3.4). Its MR constant Z remains 
undetermined, but this is not already the problem, since the renormalized gluon propagator (3.7) depends on the 
renormalized mass gap A^. So we consider it as the physical mass gap within our approach. Precisely this quantity 
should be positive, finite, gauge- and a-independent, etc., it should exist when A — > oo and e — > + . The renormaliza- 
tion of the mass gap is an example of the NP MR program. We were able to accumulate all the quadratic divergences 
into the renormalization of the mass gap. Due to the renormalization the quadratic divergences parameterized as the 
regularized mass gap A 2 may be absorbed in its re-definition leading thus finally to the physical mass gap. 

Concluding, since all the parameters in the expansion (3.7) are expressed in the IR renormalized quantities, we can 
go to the e — > + limit, without encountering any problems in this limit now. Also, we have found of no practical use 
to introduce the MR constant Z 3 for the gluon propagator itself separately from that of the mass gap (3.4). 

A. The general criterion of gluon confinement 

We are now in the position to analytically formulate the general criterion of gluon confinement. It will be instructive 
to analytically formulate the general criterion of quark confinement as well. Let us begin with the former one. 

The dimensionally renormalized (within the DT complemented by the DRM) expression for the relevant gluon 
propagator in INP QCD is Eq. (3.7). Substituting back into this equation the expansion (3.2), one obtains 



D^(g) = e x ^T^) ]T(AiO fe+1 ^ 
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(3.9) 
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Let us underline that this is the general expression for the renormalized gluon propagator, since it does not depend 
on whether the gluon momentum is independent skeleton loop variable or not. Due to the distribution nature of the 
NP IR singularities (q 2 )~ 2 ~ k , k — 0,1,2,3.... (see appendix A as well), the two principally different cases should be 
separately considered. 

1. If the gluon momentum q is an independent skeleton loop variable, then, as emphasized repeatedly above, the 
initial (q 2 )~ 2 ~~ k NP IR singularities should be regularized with the help of the expansion (3.2). Finally one arrives at 
Eq. (3.7) in the e — > + limit, as it should be. 

Let us note in advance that beyond the one-loop skeleton integrals the analysis should be done in a more sophisti- 
cated way, otherwise the appearance of the product of at least two S functions at the same point is possible. However 
this product is not defined in the DT Q ■ So in the multi-loop skeleton diagrams instead of the S functions in the 
residues their derivatives may appear (seeRef. 0] and appendix A in this paper). They should be treated in the sense 
of the DT. Fortunately, as mentioned in appendix A, the IR rcnormalization of the theory is not undermined, since 
a pole in e is always a simple pole 1/e for each independent skeleton loop variable (see the dimensionally regularized 
expansion (A10)). That is why the starting expression for the relevant gluon propagator in INP QCD is always the 
general expression (3.9). It makes it possible to perform formal algebraic operations (first of all multiplication) on 
the corresponding Laurent expansions. The final product is again the Laurent expansion. This allows one to use the 
dimensionally regularized expansion (A10), and thus to avoid the multiplication of the (5-functions at the same point, 
if the starting expansion will be expansion (3.7) for each gluon propagator. 

2. The necessary and sufficient conditions for gluon confinement are: 

(i). If, however, the gluon momentum q is not a loop variable (i.e., it is external momentum), then the initial 
{q 2 )~ 2 ~ k NP IR singularities cannot be treated as distributions, i.e., the regularization expansion (3.2) is not the case 
to be used. The functions (q 2 )~ 2 ~ k are the standard ones, and the relevant gluon propagator (3.9) vanishes as e goes 
to zero, i.e., 



oc 

D* (g) = e x iT^q) ^(A 2 ^ V)~ 2 ~ fc ~ C ^ 0+. (3.10) 

It is worth emphasizing that the final e — ► + limit is permitted to take only after expressing all the Green's functions, 
parameters and mass gap in terms of their IR renormalized counterparts because they, by definition, exist in this 
limit. This behavior is gauge-invariant, does not depend on any truncations/approximations, etc., and thus it is a 
general one. It prevents the transversal dressed gluons to appear in asymptotic states, so color dressed gluons can 
never be isolated. This is the first necessary condition for gluon confinement. 

(ii). The second sufficient condition of gluon confinement is the absence of the free gluons in the corresponding 
theory. Just such a theory has been formulated in our previous works 1, 2], namely INP QCD. Let us remind that 
its full gluon propagator (2.1) has no free gluon propagator limit when the interaction is to be switched off (A 2 = 0). 
We argued that it should be used for the calculations of physical observables, processes in low-energy QCD from first 
principles. We consider both the suppression of the colored dressed gluons at large distances and the absence of the 
free colored gluons in the theory as the exact criterion of gluon confinement (for its initial formulation see Ref. Q). 



B. The general criterion of quark confinement 

It is instructive to formulate here the quark confinement criterion in advance as well. It consists of the two 
independent conditions. 

1. The first necessary condition, formulated at the fundamental (microscopic) quark-gluon level, is the absence 
of the pole-type singularities in the quark Green's function at any gauge on the real axe at some finite point in the 
complex momentum plane, i.e., 



Sip) . (3.11) 
p- m ph 

where is the standard quark wave function renormalization constant, while m p h is the mass to which a physical 
meaning could be assigned. In other words, the quark always remains an off-mass-shell object. Such an understanding 
(interpretation) of quark confinement comes apparently from the Gribov's approach to quark confinement II 01 and 
Preparata's massive quark model (MQM) in which external quark legs were approximated by entire functions A 
quark propagator may or may not be an entire function, but in any case the pole of the first order (like the electron 
propagator has in QED) should disappear (see, for example Refs. [II, [01, [lj| an d references therein). 



5 



2. The second sufficient condition, formulated at the hadronic (macroscopic) level, is the existence of the 
discrete spectrum only (no continuum) in bound-states, in order to prevent quarks to appear in asymptotic states. 
This condition comes apparently from the 't Hooft's model for two-dimensional QCD with large N c limit [l5[ (see 
alsoRefs. pi El). 

This definition of quark confinement in the momentum space is gauge invariant, flavor independent, i.e., valid for 
all types of quarks (light or heavy), etc., and thus it is a general one. The Wilson criterion of quark confinement 
formulated in the configuration space - Area law [I?] . [l8| is relevant only for heavy quarks, as well as a linear rising 
potential between static (heavy) quarks pj|, also "seen" by lattice QCD [20I l2l"l| . 

At nonzero temperature and density, for example in the quark-gluon plasma (QGP) [H[23| (and references therein), 
the bound-states will be dissolved, so the second sufficient condition does not work any more. However, the first 
necessary condition remains always valid, of course. In other words, by increasing temperature or density there is no 
way to put quarks on the mass-shell. So that what is known as the De-confinement phase transition in QGP is in fact 
the De-hadronization phase transition. De-confinement is about the liberation of the colored objects from the vacuum 
and not from the bound-states. In the QCD ground-state there are many colored objects such as quarks, gluons, 
instantons and may be something else. Since color confinement is absolute and permanent, none of these colored 
objects can appear in physical spectrum, and thus De-confinement phase transition does not exist, in principle. 

IV. PHYSICAL LIMITS 

We introduced the renormalized mass gap in the relation (3.4), defined its existence when the dimensionless UV 
regulating parameter A goes to infinity, i.e., in the A — > oo limit. However, nothing was said about the behavior of 
the coupling constant squared g 2 in this limit. In principal, it may also depend on A, becoming thus the so-called 
"running" effective charge g 2 ~ a s = a s (X). In the general composition (3.4) 

Z- 1 (X,a s (X))A 2 (X,a s (X)), (4.1) 

all the possible types of the effective charge behavior in the A — > oo limit should be considered independently from 
each other (the dependence on other parameters is not shown explicitly, as unimportant for further discussion). 

1) . If a s (X) — > oo as A — > oo, then one recovers the strong coupling regime. Evidently, just this finite limit can be 
defined as the renormalized mass gap (3.4), i.e., in fact 

Z- 1 (A,a s (A))A 2 (A,a s (A)) = A 2 R , X -> oo, a s (X) -> oo. (4.2) 

App arently, only this mass gap can be identified/related with/to the Jaffe - Witten (JW) mass gap discussed in Ref. 
|24j (see section VII below as well). 

2) . If a s (X) — > c as A — > oo, where c is a finite constant, then it can be put unity, not losing generality. This means 
that the effective charge becomes unity, and this is only possible for the free gluon propagator. But the free gluon 
propagator contains none of the mass scale parameters, so in fact 

Z- 1 (A,a s (A))A 2 (A,a s (A)) = 0, A -> oo, a s (X) -> 1. (4.3) 

3) . If a s (A) — > as A — ► oo, then one recovers the weak coupling regime. Evidently, just this finite limit can be 
defined as Ag C£) , (however, see section VII below) i.e., in fact 

Z- 1 (A,a s (A))A 2 (A,a s (A)) = A 2 QCD , X - oo, a. (A) -> 0. (4.4) 

There is no doubt that the regularized mass gap may provide the existence of the two different physical mass scale 
parameters after the renormalization program is performed. Though these two physical parameters show up explicitly 
at different regimes, nevertheless, numerically they may not be very different, indeed, as emphasized in Ref. 
Our mass gap A 2 ^ determines the power-type deviation of the full gluon propagator from the free one in the q 2 -> 
limit. The region of small q 2 is interesting for all the NP effects in QCD. This once more emphasizes the close link 
between the behavior of QCD at large distances and its INP dynamics. At the same time, the asymptotic QCD scale 
parameter Aq CD determines much more weaker logarithmic deviation of the full gluon propagator from the free one 
in the q 2 — > oo limit (see subsection A below and Refs. 0,1, i,i). Then an interesting question arises within our 
approach. How does exactly the regularized mass gap provide the appearance of Aq CD under the PT logarithms? The 
problem is that in the full gluon propagator the regularized mass gap contribution is linearly suppressed in comparison 
with the logarithmical divergent term in the PT q 2 — > oo limit (see Refs. [H, 0). 
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A. Asymptotic freedom 

In the PT effective charge d PT (q 2 ) = [l+II s (g 2 ; d PT ))~ 1 the invariant function H s (q 2 ; d PT ) can be only logarithmical 
divergent [J, Q in the PT q 2 — > oo limit at any d, in particular at d = d PT . So putting for further convenience 
d PT (q 2 ) — a s (q 2 ; A 2 )/a s (X) in this relation, one obtains 

and b is the standard color group factor. This expression represents the summation of the so-called main PT logarithms 
in powers of a s (A). However, nothing should depend on A (and hence on A) when they go to infinity in order to 
recover the finite effective charge in this limit. To show explicitly that this finite limit exists, let us formally write 

A 2 = /(A)A 2 (A,a s (A)), (4.6) 

which is always valid, since /(A) is, in general, an arbitrary dimensionless function. In this connection, let us again 
remind that in order to get the expression (4.5) from the full effective charge d(q 2 ) = [1 + H s (q 2 ; d) + c(d)(A 2 /g 2 )] -1 
0, IH the mass gap contribution A 2 /q 2 is only asymptotically suppressed in the q 2 — > oo limit. In other words, we 
distinguish between the asymptotic suppression of the mass gap contribution A 2 jq 2 and the formal PT A 2 = limit. 
So the mass gap A 2 = A 2 (A, a s (X)) itself here is not put identically zero and hence the relation (4.6) makes sense. 
On account of the relation (4.4), it becomes 

A 2 = f(X)Z(X,a s (X))A 2 QCDl A^co, a s (X) -> 0. (4.7) 
Substituting it into the expression (4.5) and doing some algebra, one obtains 



1 + ba s \n(q 2 /A z QCD ) 



if and only if 

a. (A) 



1 - ba s (X) In(fZ) 



, A^oo, a s (A)->0. (4.9) 



exists and is finite in the above shown limits. Here we introduce the short hand notations / = /(A) and Z = 
Z(X,a s (X)), for simplicity. Evidently, the finite a s can be identified with the fine-structure constant of the strong 
interactions, calculated at some fixed scale. It is worth emphasizing that the existence and finiteness of a s is due to 
the product (fZ). Indeed, from Eq. (4.9) it follows 

Hf z ) = ^Stttt ~> r^Tiv A ^°°< «*( A )-°> ( 4 - 10 ) 

a s ba s {X) ba s (X) 

which means that (fZ) — exp(l/6a s (A) in the above shown limits. Substituting this into the relation (4.7) it becomes 

(Ajfcoo A2 ° XP (-^(A)) - ^QCDi «.(A) - 0, (4.11) 

which is the finite limit of the renormalization group equations solution [25[ . 
At very large q 2 from Eq. (4.8) one recovers 



which is nothing but asymptotic freedom (AF) famous formula if b > @, SL H & E2 ■ In QCD with three colors and 
sixth flavors this is so, indeed. For the pure Yang-Mills (YM) fields b = 11 /An > always. Let us underline that in 
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the expressions (4.8) and (4.12) q 2 is always big enough, so it cannot go below Aq CD . We have shown explicitly the 
AF behavior of QCD at short distances (q 2 — ► oo), not using the renormalization group equations and their solutions 
(i.e., we need no expansion in powers of the coupling constant for the corresponding /3-function) 0, 0, H, B HH ■ The 
regularized mass gap is suppressed in the q 2 — > oo limit, as it has been mentioned above. From Eq. (4.4) and our 
consideration in this subsection it follows, nevertheless, that the regularized mass gap in the A — > oo limit provides 
the existence of the asymptotic QCD scale parameter Aq C£) as well. 

There is no relation between the renormalized mass gap A 2 R (4.2) and the asymptotic scale parameter ^q CD (4.4), 
since they show up explicitly at different regimes. They are different scales, indeed, responsible for different NP and 
non-trivial PT dynamics in QCD, though numerically they may not be very different, as underlined above. However, 
originally they have been generated in the region of small q 2 . In Ref. (2o| it has been noticed that being numerically 
a few hundred MeV only, Aq CD cannot survive in the q 2 — > oo limit. So none of the finite mass scale parameters can 
be determined by the PT QCD. They should come from the region of small q 2 , being thus NP by origin and surviving 
the renormalizatin program (i.e, the removal of A in the A — > oo limit), as was just demonstrated above. 

Concluding, all this can be a manifestation that "the problems encountered in perturbation theory are not mere 
mathematical artifacts but rather signify deep properties of the full theory" [26(. The message that we are trying to 
convey is that the INP dynamical structure of the full gluon propagator indicates the existence of its nontrivial PT 
one and the other way around. 



V. CONFINING POTENTIAL 



After discussing some important aspects of the general MR program for the INP gluon propagator using the DRM 
and the DT, it makes sense to go back to the initial expansion for the INP gluon propagator (2.1). A new surprising 
feature of this solution is that its structure at zero (q 2 — > 0) can be again determined by its A 2 /(q 2 ) 2 term only. To 
show this explicitly, let us begin with the theorem from the theory of functions of complex variable [27| , which is an 
extremely useful for the explanation of the behavior of our solution (2.1) in the deep IR q 2 — > limit. 

The function L(q 2 ) is defined by its Laurent expansion (2.2), and thus it has an isolated essentially singular point 
at q 2 = 0. Its behavior in the neighborhood of this point is regulated by the Weierstrass-Sokhatsky-Casorati (WSC) 
theorem (see appendix B) which tells that 

]im L(q 2 n )=Y, q 2 n ^0, (5.1) 

n — >oo 

where Y is any complex number, and q 2 is a sequence of points q\ , g|...q 2 ... along which q 2 goes to zero, and for which 
the above-displayed limit always exists. Of course, Y remains arbitrary (it depends on the chosen sequence of points 
q 2 ), but, in general, it depends on the same set of parameters as the residues, i.e., Y = Y(\,a,t;,g 2 ). This theorem 
thus allows one to replace the Laurent expansion L(q 2 ) by Y when q 2 — > independently from all other test functions 
in the corresponding skeleton loop integrands, i.e., 

L{q 2 )^Y(\a^ l9 2 ), q 2 ^0. (5.2) 

Our consideration in this section up to this point is necessarily formal, since the mass gap remains unrenormalized. 
So far it has been only regularized, i.e., A 2 = A 2 (A, a, £, g 2 ). The renormalization of the mass gap can be proceed as 
follows. Due to the above-formulated WSC theorem, the full gluon propagator (2.1) effectively becomes 

D^(q) = D^q;A 2 )=iT^q)-^Y(\a,^g 2 )A 2 (\a^,g 2 ), q 2 -> 0, (5.3) 

so just the A 2 (g 2 ) -2 structure of the full gluon propagator (2.1) is all that matters, indeed. Let us now define the 
renormalized (R) mass gap as follows: 



A R = Y(\,a,Z,g 2 )A 2 (\,a,Z,g 2 ), (5.4) 

so that we consider Y(X, a, £, g 2 ) as the MR constant for the mass gap, and A 2 R is the physical mass gap within our 
approach. Precisely this quantity should be positive, finite, gauge- independent, etc., it should exist when A — > oo and 
a — > O.Due to the WSC theorem, we can always choose such Y = Y\ x Y" 2 x ••■ in order to satisfy all the necessary 
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requirements (each Y n will depend on its own sequence of points along which q 2 — ► 0, the so-called subsequences). 
Numerically it should be identified with the mass gap A 2 R introduced by the relation (3.4) and described in section 
III before subsection A (see also the relation (5.8) below). 

Thus the full gluon propagator relevant for the INP QCD dynamics becomes 



D^(q;A 2 R )=iT IM/ (q)-^. (5.5) 

It is possible to say that because of the WSC theorem there always exists such a sequence of points along which q 2 — > 
that the Laurent expansion (2.1)-(2.2) effectively converges to the function (5.5) in the whole <7 2 -momentum plane 
(containing both points q 2 = and q 2 = oo) after the renormalization of the mass gap is performed. But the region 
of small q 2 is of special interest. In this region the confinement dynamics begin to play a dominant role. However, 
it is worth emphasizing that the solution (5.5) is not the IR asymptotic (q 2 —> 0) of the initial Laurent expansion 
(2.1)-(2.2), since {q 2 )~ 2 term is not a leading one in this limit, it is a leading one just in the opposite q 2 — > oo limit. 

The severe IR singularity (q 2 )~ 2 , which only one is present in the relevant gluon propagator (5.5), is the first NP IR 
singularity possible in four-dimensional QCD. A special regularization expansion is to be used in order to deal with 
it, since the standard methods fail in this case. As mentioned above, it should be correctly treated within the DT Q 
complemented by the DRM Q. If q 2 is an independent skeleton loop variable, then the corresponding dimensional 
regularization of this singularity is given by the expansion (3.2) at k = 0, namely 



(<7 2 )- 



l r 

e 



a(0)S 4 (q) + O a (e) 



H 



(5.6) 



and a(0) = tt 2 . Due to the S 4: (q) function in the residue of this expansion, all the test functions which appear under 
corresponding skeleton loop integrals should be finally replaced by their expression at q = 0. So the dimcnsionally 
regularized expansion for the gluon propagator (5.5) becomes 



D^(q; A|) = iT^(q)^R * 1 a(0)SHq) + 0(e) 



0" 



(5.7) 



where we put Oo(e) = O(e), for simplicity. As emphasized above, in the presence of the IR regularization parameter 
e all the Green's functions and parameters become, in general, dependent on it. The only way to remove the pole 1/e 
from the relevant gluon propagator (5.7) is to define the IR renormalized mass gap as follows: 



A 2 R = X(e)A 2 R = exA 2 R , e^0+, (5.8) 

where X(e) — e is the IR MR constant for the mass gap. Let us remind that contrary to A 2 R its IR renormalized 
counterpart A^, exist as e — > + . In both expressions for the mass gap the dependence on e is assumed but not shown 
explicitly. So we distinguish between both mass gaps only by the dependence on e. 

On account of the relation (5.8), the dimensionally regularized expansion (5.7) finally becomes 



D^(q; A 2 R ) = tT^(q)A 2 R a(0)6 4 (q) + 0(e), e 0+ (5.9) 

Evidently, after performing the renormalization program (i.e., going to the IR renormalized quantities), the terms of 
the order 0(e) can be omitted from the consideration. 

The renormalization of the mass gap automatically IR renormalizes the relevant gluon propagator (5.5), so that 



D^(q-,A 2 R )=exiT^(q)j-^, e - 0+, (5.10) 

in complete agreement with the expansion (3.9) at k = with putting there CIq = 1. In order to achieve the agreement 
with the general renormalization program performed in section III it is necessary to put Y~ (A, e)X(e) = Z(X, e) and 
hence $o(A; e) — CIq Y(\, e) = Y(X, e), on account of the convergence condition (3.6) at k = and the relation (5.8). 

There is no doubt left that the WSC theorem somehow underlines the importance of the simplest NP IR singularity 
l/(q 2 ) 2 possible in four-dimensional QCD, while all other may be suppressed in the deep IR (q 2 — > 0) region due 
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to this theorem. However, from the general consideration in sections III and IV the suppression mechanism is not 
seen, at least at this stage. In principle, one can develop the formal PT series in powers of the mass gap. Then the 
first nontrivial approximation is just the above-mentioned singularity A R /(q 2 ) 2 , which is only one to be started with, 
anyway. It is well known that the expression (5.5) leads to the linear rising potential between heavy quarks also 
"seen" by lattice QCD [2(J[2lj]. That is why we call it the confining potential. Just it will be used for the derivation 
of the system of equations determining the confining quark propagator (for preliminary consideration see Ref. [28| 
and references therein). Evidently, what we have pointed out in subsection A of section III for the general case is 
valid for the confining potential (5.5) as well. 



A. The renormalized "running" effective charge 

It is instructive to find explicitly the corresponding /3-function. From Eq. (5.5) it follows that the corresponding 
Lorentz structure, which is nothing but the corresponding effective charge ("running"), in terms of the renormalized 
mass gap is 



d(q 2 ;Al) = a s (c?;Al) = ^k, (5.11) 

and this does not depend whether the gluon momentum is independent loop variable or not. Then from the renor- 
malization group equation for the renormalized effective charge, which determines the /3-function, 



dq 

it simply follows that 



f3(a s (q 2 ; A 2 ,)) = -a s (q 2 ; A 2 B ) = --f (5.13) 

Thus, one can conclude that the corresponding /3-function as a function of its argument is always in the domain of 
attraction (i.e., negative). So it has no IR stable fixed point indeed as it is required for the confining theory [3J. Just 
these expressions for the /3-function and the running effective charge should be used for the calculation of the truly 
NP quantities in the YM theory, such as the gluon condensate, the gluon part of the Bag constant p9l. [30l[3l|. etc. 
In phenomenology for these purposes we need to know the ratio j3(q 2 )/a s (q 2 ) in the q 2 — > limit, so it is always 



P(g 2 ) _ P(a s (q 2 ;A 2 R )) 

a s (q 2 )- a s (q 2 ;A 2 R ) > 

within our approach to low-energy QCD, that's INP QCD (don't forget that the replacement a s (q 2 ) — > al NP (q 2 ) and 
f3(q 2 ) -» f3 INP (q 2 ) is assumed @). 

Concluding, let us note if the renormalized effective charge (5.11) is to be multiplied by the additional powers of 
(q 2 )~ 2 ~ k , k — 0, 1,2,3, then this product should be treated as the full gluon propagator itself, described above. 
The only difference between them becomes the unimportant tensor structure T^^q). 



VI. DISCUSSION 



Thus the QCD vacuum is really beset with severe IR singularities, which have been summarized (accumulated) 
into the INP gluon propagator (2.1). There is no doubt that the purely transversal severely singular virtual gluon 
field configurations play important role in the dynamical and topological structure of the QCD ground state, leading 
thus to the general zero momentum modes enhancement (ZMME) effect there reflected in the INP gluon propagator. 
Evidently, the ZMME (or simply ZME) mechanism of confinement (see our previous papers (321 . |33| and references 
therein as well) is nothin g bu t the well forgotten IR slavery (IRS) one, which can be equivalcntly referred to as a 
strong coupling regime [S Il6|. 
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Indeed, at the very beginning of QCD it was expressed a general idea [H, [H, HH, [H, H3, [H, H^, H(| that the 
quantum excitations of the IR degrees of freedom, because of self-interaction of massless gluons in the QCD vacuum, 
made it only possible to understand confinement, dynamical (spontaneous) breakdown of chiral symmetry and other 
NP effects. In other words, the importance of the deep IR structure of the QCD vacuum has been emphasized as well 
as its relevance to the above-mentioned NP effects and the other way around. This development was stopped by the 
wide-spread wrong opinion that severe IR singularities cannot be put under control. Here we have explicitly shown 
(see also our recent papers [H, l4lj and references therein) that the adequate mathematical theory for quantum YM 
gauge theory is the DT (the theory of generalized functions) [7[, complemented by the DRM [8|. Together with the 
theory of functions of complex variable [13] they provide a correct treatment of these severe IR singularities without 
any problems. Thus, we come back to the old idea but on a new basis that is why it becomes new ("new is well 
forgotten old"). In other words, we put the IRS mechanism of color confinement on a firm mathematical ground. 
This makes it possible to analytically formulate the gluon confinement criterion in a gauge invariant way for the first 
time. 

The confining potential (5.5) in the different approximations and gauges has been earlier obtained and investigated 
in many papers (see, for example Ref. [42 | and references therein). We have confirmed and thus revitalized these 
investigations, in which this behavior has been obtained as an IR asymptotic solution to the gluon SD equation. 
However, let us emphasize once more that due to the WSC theorem the confining potential (5.5) is not the IR 
asymptotic of the initial Laurent expansion (2.1)-(2.2). Moreover, the whole INP gluon propagator (2.1) effectively 
converges to the confining potential (5.5) after the renormalization of the mass gap is performed. The WSC theorem 
clearly shows that the Aj i /(q 2 ) 2 structure is only important, while all other terms in the INP gluon propagator (2.1) 
are suppressed (though each next term in the expansion (2.2) is more singular in the IR than the previous one). 

In the presence of the mass gap the cou plin g constant plays no role. This is also an evidence of the "dimensional 
transmutation" , g 2 — > A 2 (A, a, £, g 2 ) [H. l43l. |44| . which occurs whenever a massless theory acquires masses dynamically. 
It is a general feature of spontaneous symmetry breaking in field theories. In our case, the color gauge symmetry is 
broken at the level of the full gluon self-energy, while maintaining at the level of the full gluon propagator, which is more 
important [l[ (something like a "self-consistent violation" or a "hidden violation" of the color gauge symmetry). In the 
massive solution |2j the mass gap transforms further into the effective gluon mass, i.e., g 2 — > A 2 (A,a,£, g 2 ) — ► m ff (£). 
Nevertheless, it remains gauge-dependent (i.e., not physical) even after the corresponding full gluon propagator 
is renormalized. Within the INP QCD the mass gap transforms further into the physical mass gap, i.e., g 2 — > 
A 2 (A, a, £, g 2 ) — > Ap, and the gluons remain massless in a gauge invariant way (this paper). Let us note that the 
renormalization of the mass gap automatically renormalizes the INP QCD gluon propagator as well. 

VII. CONCLUSIONS 

Let us denote the version of our mass gap Ap which will appear in the S-matrix elements for the corresponding 
physical quantities/processes in low-energy QCD as A 2 NP (in principle, they may be slightly different from each other, 
indeed). Then a symbolic relation between it and the initial mass gap A 2 (A,a s (A)) and Ap T instead of Aq CD (for 
reason see discussion below) could be written as follows: 

ALvp ^H s(A) A 2 (A,a s (A)) — A 2 PT , (7.1) 

which summarizes our main results in this investigation. QCD as a quantum gauge field theory, describing the 
interactions of never seen colored objects (the gluons and quarks), cannot have the physical mass gap. In other words, 
this is a theory which describes the behavior of the colored objects in the vacuum. In QCD mass gap may only 
appear in the way described in our previous publication [l[, that's A 2 (A, a s (X)) in Eq. (7.1). In order to become the 
theory of the strong interactions it should undergo the two phase transitions; in the weak and strong coupling regimes. 
In the first case it becomes the PT QCD which describes all the high-energy phenomena in the strong interactions 
from first principles (AF, scale violation, hard processes, etc.). It has its own physical mass gap which we denote as 
A|, T in Eq. (7.1). In the second case it becomes the INP QCD which describes all the low-energy phenomena in the 
strong interactions from first principles (those includes first of all color confinement, dynamical breakdown of chiral 
symmetry, bound-states, etc.). It has its own physical mass gap, that's A 2 INP in Eq. (7.1). 

In this connection, a few things should be made perfectly clear. First of all, let us underline that the PT QCD and 
the INP QCD are not effective theories, as pointed out above both theories are fundamental ones. Secondly, such a 
quantity as Aq CD does not exist at all, since QCD itself cannot have a physical limit. In order to avoid any confusion 
the corresponding scale is better to denote as Ap T , since just the PT QCD is responsible for all the high-energy 
phenomena in the strong interactions. Thus similarly the relation (7.1), the following symbolic relation makes sense 
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INP QCD QCD PT QCD, (7.2) 

so that at the fundamental (quark-gluon) level the PT QCD is AF, while the INP QCD confines gluons. In the 
subsequent papers we will show that this theory will confine quarks, as well as will explain the spontaneous breakdown 
of chiral symmetry. Both theories have their own mass gaps A 2 wf , and Ap T , which are solely responsible for the large- 
and short-scale structures of the QCD ground state, respectively. 

A few years ago Jaffe and Witten have formulated the following theorem [2l[ : 

Yang-Mills Existence And Mass Gap: Prove that for any compact simple gauge group G, quantum Yang-Mills 
theory on R 4 exists and has a mass gap A > 0. 

Of course, to prove the existence of the YM theory with compact simple gauge group G is a formidable task yet. 
It is rather a mathematical than a physical problem. However, from the JW presentation of their theorem it clearly 
follows that their mass gap should be identified with our mass gap K 2 INP . At the same time, we have argued above that 
QCD itself cannot have a physical mass gap. It has a mass gap which is only regularized, i.e., A 2 = A 2 (A, a s (X)), and 
therefore there is no guarantee that it is positive. It cannot be related directly to any of physical quantities/processes. 
Let us also remind that QCD cannot confine free gluons |l|, |2| . As actual theory of the strong interactions the two 
different faces of QCD come into the play: the PT QCD for high-energy physics and INP QCD for low-energy physics, 
which confines "dressed" gluons, while free gluons do not exist in this theory. The corresponding mass gaps have now 
physical meanings (they are finite, positive, gauge-invariant, etc.). 

Our basic result obtained in the previous works [l|, 0] and in this paper can be jointly formulated as follows: 

Mass Gap Existence And Gluon Confinement: If quantum Yang-Mills theory with compact simple gauge 
group G = SU(3) exists on R 4 , then undergoing the phase transition in the strong coupling regime it becomes INP 
QCD, which has a physical mass gap and confines gluons. 

Some important features of the INP QCD are: 

1. Its full gluon propagator (2.1) converges to the expression (5.5) after the renormalization of the mass gap is 
performed. This expression is effectively valid in the whole q 2 -momentum plane. 

2. It has a physical mass gap. 

3. It confines "dressed" gluons in asymptotic states. 

4. It has no free gluons at all due to the subtraction method proposed and formulated in our previous works [3, [2]. 
Some other interesting features of this theory may be established after the explicit including of the quark degrees 

of freedom into its formalism in the next our papers. However, first in the subsequent paper we will show how the 
QCD vacuum structure is to be investigated within the YM version of the INP QCD. 

Concluding, a few remarks are in order. In our previous work [l|, we have explicitly shown that the full photon 
propagator in quantum electrodynamics (QED) has only the PT-type IR singularity, 1/g 2 . This is in agreement with 
the cluster property of the Wightman functions |45j, that's correlation functions of observables. In QCD the explicit 
presence of the regularized mass gap, which are necessarily accompanied by severe IR singularities (g 2 ) -2 ~ fe , k = 
0,1,2,3..., apparently, will violate this property. In turn, this validates the Strocchi theorem [46], which allows 
such a severely singular behavior of the full gluon propagator in QCD. However, this is not a problem, since QCD 
has no physical observables. In PT QCD, which gluon propagator is as much singular as l/q 2 only, the cluster 
property will not be violated. On the other hand, in INP QCD with a such singular behavior of the relevant gluon 
propagator (5.5) the situation with the Wightman functions is not clear. It can be clarified only after the solution of 
the color confinement problem, and a realistic calculations of the various physical observables within INP QCD. At 
the fundamental quark-gluon level only these remarks make sense about the correlation between the structure of the 
corresponding gluon propagator and the properties of the Wightman functions. 
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APPENDIX A: INFRARED DIMENSIONAL REGULARIZATION WITHIN THE DISTRIBUTION 

THEORY 



1. The DRM in the PT 



As repeatedly emphasized in our previous works [l|, [2| , the mass gap in nothing but the re-defined skeleton tadpole 
term, or it can be reduced to the skeleton tadpole term itself. Also, it has been explicitly shown that there is no such 
regularization scheme (preserving or not gauge invariance) in which the transversality condition for the full gluon 
self-energy could be satisfied unless the constant skeleton tadpole term 



U^D) = U t (D) ee A 2 t (D) = g 2 J ^T°D( gi ), (Al) 

is to be disregarded from the very beginning, i.e., put formally zero everywhere. Here T4 is the four-gluon point-like 
vertex, and g 2 is the dimensionless coupling constant squared. We omit the tensor and color indices in this integral, 
as unimportant for further discussion. It is nothing else but the quadratically divergent in the PT constant, so it is 
assumed to be regularized. The mass gap is not survived in the PT q 2 — > 00 limit |l|, however, in the PT there are 
still problems with such kind of integrals. 

In the PT, when the first non-trivial approximation for the full gluon propagator D = D(q) is the free one Do = 
Do(q), the constant tadpole term is to be simply discarded, i.e., to be put formally zero within the DRM 0], so 
that II^Do) = S p a A% (Do) = 0. However, even in the DRM this is not an exact result, but rather an embarrassing 
prescription, as pointed out in Ref. ||. To show explicitly that there are still problems, as mentioned above, it is 
instructive to substitute the first iteration of the gluon SD equation into the previous expression (Al). Symbolically 
it looks like D(q) = D (q) + D (q)iU(q; D)D(q) = D (q) + D (q)iU(q; D )D (q) + ... = D (q) + D^(q) + where 
we omit all the indices and put Do = D(°\ Doing so, one obtains 



U t (D = D + D^ + ...) = U t (Do)+H t (D^) + ...^U t (D )+g 2 J l A^T°[D ( qi )] 2 m( qi ; D ) + ... 

= n t (D )+Tl t (Do)g 2 J ^T°[Do( qi )] 2 +g 2 J ^T°[Do( qi )?W (q x] Do) + ... . 

(A2) 

Here we introduce the subtraction as follows: H s (qi;Do) = H(qi;Do) — H(0;Do), and put H(0; Do) = lit (Do), for 
simplicity, when the mass gap is to be reduced to the tadpole term itself [l|. In the second line of Eq. (A2) the 
first integral is not only UV divergent but IR singular as well. If we now omit the first term in accordance with 
the above-mentioned prescription, the product of this integral and the tadpole term H t (Do) remains, nevertheless, 
undetermined. Moreover, the structure of the second integral in this line is much more complicated than in the 
divergent constant integral n t (Do) in Eq. (Al). All this reflects the general problem that such kind of massless 
integrals 



d a q gMi-g/ip , A on 
(2ir) d (q 2 ) n 1 > 

are ill defined, since there is no dimension where they are meaningful. They are either IR singular or UV divergent, 
depending on the relation between the numbers d, p and n ||. This prescription clearly shows that the DRM, 
preserving gauge invariance, nevertheless, is by itself not sufficient to provide us insights into the correct treatment 
of the power-type IR singularities shown in Eq. (A3) (we will address this problem below in subsection 2). Thus, one 
concludes that the tadpole term (Al) A 2 (D) = A 2 (X, a;D) is, in general, not zero. 

However, in the PT we can adhere to the prescription that such massless tadpole integrals can be discarded in the 
DRM 0, [§]. As mentioned above, we have already shown [l[ that the mass gap, in general, and the skeleton tadpole 
term, in particular, can be neglected in the PT, indeed (not depending on whether A, a are to be introduced within 
the regularization scheme preserving gauge invariance or not). In what follows we will show how precisely the DRM 
Q should be correctly implemented into the DT Q in order to control the power-type severe IR singularities, which 
may appear not only in the PT series, but mainly in the NP QCD as well (see Ref. Q). 
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2. The DRM in the DT 



In general, all the Green's functions in QCD are generalized functions, i.e., they are distributions. This is especially 
true for the NP IR singularities of the full gluon propagator due to the self-interaction of massless gluons in the QCD 
vacuum. They present a rather broad and important class of functions with algebraic singularities, i.e., functions with 
nonsummable singularities at isolated points [7[ (at zero in our case). Roughly speaking, this means that all relations 
involving distributions should be considered under corresponding integrals, taking into account the smoothness prop- 
erties of the corresponding space of test functions. Let us note in advance that the space in which our generalized 
functions are continuous linear functionals is K, that's the space of infinitely diffcrcntiable functions having compact 
support, i.e., they are zero outside some finite region (different for each differentiable function) Q. 

Let us consider the positively definite (P > 0) squared (quadratic) Euclidean form 



P(q) = <?o + <?i + + - + in-i = r, (A4) 

where n is the number of the components. The generalized function (distribution) P x (q), where A being, in general, 
an arbitrary complex number, is defined as 



(P A ,^) = [ P\qMq)d d q, (A5) 
Jp>a 

where tp(q) is the above-mentioned some test function. At ReX > this integral is convergent and is an analytic 
function of A. Analytical continuation to the region ReX < shows that it has a simple pole at points [7fl 



X = ---k, k = 0,1,2,3... (A6) 

In order to actually define the skeleton loop integrals in the deep IR domain, which appear in the system of the 
SD equations, it is necessary to introduce the IR regularization parameter e, defined as d = n + 2e, e — ► + within 
the DRM [H, where d is the dimension of the loop integral (see Eq. (A5)). As a result, all the Green's functions and 
"bare" parameters should be regularized with respect to e which should be set to zero at the end of the computations. 
The structure of the NP IR singularities is then determined (when n is even number) as follows [3] : 



c (k) 

^ X = xTwWTk +fimteterms > (A7) 



where the residue is 



n /2 

^ = 2»fc!r«n/2) + fc) X LkSn{q) (A8) 

with L = (d 2 /dql) + (dydqf) + ... + (P/dq*^). 

Thus the regularization of the NP IR singularities (A5), on account of (A6), is nothing but the whole expansion in 
the corresponding powers of e and not the separate term(s). Let us underline its most remarkable feature. The order 
of singularity does not depend on A, n and k. In terms of the IR regularization parameter e it is always a simple pole 
1/e. This means that all power terms in Eq. (A7) will have the same singularity, i.e., 



( g 2 )-f- fc = -CW + f inite terms, e -► 0+, (A9) 

where we can put d = n now (i.e., after introducing this expansion). By "finite terms" here and everywhere a 
number of necessary subtractions under corresponding integrals is understood [3]. However, the residue at a pole will 
be drastically changed from one power singularity to another. This means different solutions to the whole system of 
the SD equations for different set of numbers A and k. Different solutions mean, in their turn, different vacua. In 
this picture different vacua are to be labeled by the two independent numbers: the exponent A and k. At a given 
number of d(= n) the exponent A is always negative being integer if d{= n) is an even number or fractional if d{= n) 
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is an odd number. The number k is always integer and positive and precisely it determines the corresponding residue 
at a simple pole, see Eq. (A9). It would not be surprising if these numbers were somehow related to the nontrivial 
topology of the QCD vacuum in any dimensions. 

It is worth emphasizing that the structure of severe IR singularities in Euclidean space is much simpler than in 
Minkowski space, where kinematical (unphysical) singularities due to the light cone also exist [H, 0, H3] (in this 
connection let us remind that in Euclidean metrics q 2 = implies q t = and vice-versa, while in Minkowski metrics 
this is not so). In this case it is rather difficult to untangle them correctly from the dynamical singularities, the 
only ones which are important for the calculation of any physical observable. Also, the consideration is much more 
complicated in the configuration space That is why we always prefer to work in the momentum space (where 
propagators do not depend explicitly on the number of dimensions) with Euclidean signature. We also prefer to work 
in the covariant gauges in order to avoid peculiarities of the non-covariant gauges jj, |48| . |4^ |. for example, how to 
untangle the gauge pole from the dynamical one. 

In principle, none of the regularization schemes (how to introduce the IR regularization parameter in order to 
parameterize the NP IR divergences and thus to put them under control) should be introduced by hand. First of all, 
it should be well defined. Secondly, it should be compatible with the DT jjjj]. The DRM Q is well defined, and here 
we have shown how it should be introduced into the DT (complemented by the number of subtractions, if necessary). 
Though the so-called ±ie regularization is formally equivalent to the regularization used in our paper (see again Ref. 
0), nevertheless, it is rather inconvenient for practical use. Especially this is true for the gauge-field propagators, 
which are substantially modified due to the response of the vacuum (the ±ie prescription is designated for and is 
applicable only to the theories with the PT vacua, indeed [Tol [Hoj V Other regularization schemes are also available, 
for example, such as analytical regularization used in Ref. jl6l | or the so-called Speer's regularization (5j| . However, 
they should be compatible with the DT as emphasized above. Anyway, not the regularization is important but the 
DT itself. Just this theory provides an adequate mathematical framework for the correct treatment of all the Green's 
functions in QCD (apparently, for the first time the distribution nature of the Green's functions in quantum field 
theory has been recognized and used in Ref. [E3]). 

The regularization of the NP IR singularities in QCD is determined by the Laurent expansion (A9) at n = 4 as 
follows: 



(q 2 )- 2 - k = -a(k)[6%)}^ + f.t. = - \a(k)[5\q)}^ + O k (e)} , e -> 0+, (A10) 
e e L J 

where a(k) = 7r 2 /2 2fe fc!r(2 + k) is a finite constant depending only on k and [<5 4 (q)] ( - fe ' represents the fc's derivative 
of the 5-function (see Eqs. (A7) and (A8)). We point out that after introducing this expansion everywhere one can 
fix the number of dimensions, i.e., put d = n = A for QCD without any further problems. Indeed there will be no 
other severe IR singularities with respect to e as it goes to zero, but those explicitly shown in this expansion. Let us 
underline that, while the initial expansion (2.1) is the Laurent expansion in the inverse powers of the gluon momentum 
squared, the regularization expansion (A10) is the Laurent expansion in powers of e. This means that its regular part 
is as follows: f.t. = (q 2 )Z 2 ~ k + e(q 2 )Z 2 ~ k \nq 2 + 0(e 2 ), where for the unimportant here definition of the functional 
(q 2 )Z 2 ~ k see Ref. @- These terms, however, play no any role in the IRMR program which has been discussed in 
section III. The dimensionally regularized expansion (A10) takes place only in four-dimensional QCD with Euclidean 
signature. In other dimensions and/or Minkowski signature it is much more complicated as pointed out above. As it 
follows from this expansion any power- type NP IR singularity, including the simplest one at k = 0, scales as 1/e as it 
goes to zero. Just this plays a crucial role in the IR renormalization of the theory within our approach. Evidently, 
such kind of the dimensionally regularized expansion (A10) does not exist for the PT IR singularity, which is as much 
singular as (<; 2 ) -1 only. 

In summary, first we have emphasized the distribution nature of the NP IR singularities. Secondly, we have 
explicitly shown how the DRM should be correctly and in a gauge invariant way implemented into the DT. This 
makes it possible to put severe IR singularities under firm mathematical control. 



APPENDIX B: THE WSC THEOREM 

One of the main theorems in the theory of functions of complex variable [13] is the the above-mentioned WSC 
theorem. It describes the behavior of meromorphic functions near essential singularities. 

Theorem (Weierstrass-Sokhatsky-Casorati). If zo is an essential singularity of the function f(z), then for any 
complex number Z there exists the sequence of points z k — * zq, such that 
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lim f(z k ) = Z. (Bl) 

fc^oo 

So this theorem tells us that the behavior of the function f(z) near its essential singularity z is not determined, 
i.e, in fact, it remains arbitrary. It depends on the chosen sequence of points z k along which z goes to zero, that's 
Z = Z(z k ) (do not mixed this complex number with the MR constant of the mass gap in the relation (3.4)). 

Let us consider one classical example [13] ■ The function 



oo 

/« - el/z - E ^ ( B2 ) 



n'.z 

n=0 



has the above shown Laurent series about the essential singularity at z — 0, i.e., this Laurent expansion converges to 
the function f(z) everywhere apart from the point z = 0. At the same time, due to the WSC theorem the behavior 
of the function f(z) near z = 0, and hence of the Laurent expansion itself, depends on the chosen sequence of points 
Zk along which z — > 0. So let us proceed as follows: 
(i). If one chooses z k = 1/fc, k = 1, 2, 3, then 

lim f(z k ) = lim e fc = oo. (B3) 



fc^ 



(ii). If one chooses z k — — k = 1, 2, 3, then 



lim = lim e~ fc = 0. (B4) 

k — »oo — >OG 



(iii). If one chooses z k = l/lnA+ 2km, k = 0, 1,2,3, then 



lim f(z k ) = lim e ^ A + 2k ™ = A, (B5) 

— >oo fc — >oo 



where A is some finite constant. Many other examples can be found in Ref. [27[ and in other text-books on the theory 
of functions of complex variable. 

Concluding, let us make one important point perfectly clear. For example, the function 

v n— 

has no singularity at zero at all, while its Laurent expansion always has an essential singularity at z — 0, by definition. 
The equality (B6) means that this Laurent expansion converges to the above shown function f(z) in the ring which 
excludes zero point. Its region of convergence is 1 < \z\ < oo, while the behavior of any Laurent expansion near its 
essential singularity is always governed by the WSC theorem, in particular of the Laurent expansion shown in the 
right-hand-side of Eq. (B6). 

Another characteristic example is the function 

which is nothing but geometric series. The Laurent expansion, shown in the right-hand-side of this equation, converges 
to the function f(z) also in the region 1 < \z\ < oo, while the behavior of the Laurent expansion itself at z — > 
is again governed by the WSC theorem. The same is true for the Laurent expansion in Eq. (B2). It converges to 
the function exp(l/z) in the whole complex plane apart from the point z = 0, while near this point its behavior is 
uncertain, as it was described above. The message we are trying to convey is that all the equalities containing the 
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Laurent expansions should be treated carefully in accordance with the above-mentioned theorem (i.e., the region of 
convergence should be fixed clearly, otherwise the equality can be incorrectly understood) . 
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